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We present a general necessary and sufficient criterion for the possibility of a state transformation from 
one mixed Gaussian state to another of a bi-partite continuous-variable system with two modes. The class of 
operations that will be considered is the set of local Gaussian completely positive trace-preserving maps. 
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Imagine a physical device that is able to manipulate lo- 
cally the state of a composite quantum system by actions on 
its parts. Which state transformations could this device im- 
plement in principle, abstracting from experimental imperfec- 
tions? This question is particularly important in the field of 
quantum information theory [[!]], which concerns itself with 
the problem of whether a certain resource, e.g. an entangled 
quantum system in a known state, could be used to accom- 
plish an envisioned task. To be more specific, one asks for 
mathematical conditions that have to be met in order for a 
state transformation under natural constraints to be possible. 

Such a natural constraint is that the device can only im- 
plement local quantum operations supplemented by classi- 
cal communication (LOCC), as many applications in quan- 
tum information science involve spatially separated parties 
sharing entangled states. So far, when investigating transfor- 
mation criteria under LOCC all efforts have been devoted to 
the case where the involved quantum systems possess finite- 
dimensional Hilbert spaces as, e.g., qubit systems. The widely 
acknowledged result of Ref. [[|] relates the problem of the de- 
terministic transformation between pure states by means of 
LOCC to the mathematical theory of majorization. Based on 
this insight a series of further results have been found 
While the constraint to general LOCC is natural for low- 
dimensional systems, the situation is quite different for sys- 
tems with an infinite-dimensional Hilbert space, such as the 
modes of an electromagnetic field. The experimental opera- 
tions that are typically available are those that involve beam- 
splitters, phase shifters, and squeezers together with the abil- 
ity to prepare ancilla systems in a standard state such as the 
vacuum. The class of states that can be generated by these 
operations, and which is therefore particularly relevant from 
an experimental point of view, is the set of Gaussian states 
[Q-g]. Several properties of entangled Gaussian states are al- 
ready known. In particular, the problems of distillability and 
separability of Gaussian states have been investigated in great 
detail, and can actually be considered solved However, 
the general question of the local interconvertability between 
entangled Gaussian states has not been addressed before. 

This letter presents a first step towards finding tools for de- 
ciding whether a desired transformation of Gaussian states can 
be accomplished without the need of going through all phys- 
ical protocols, which can be an extremely tedious task. We 



will present a general necessary and sufficient criterion for the 
possibility of state transformations of a two-mode continuous- 
variable system. The class of allowed operations is the set 
of local Gaussian completely positive operations JlQ[], that is, 
those local operations that can be realized by means of local 
joint symplectic transformations on both the system and ar- 
bitrary appended ancilla systems that have been prepared in 
Gaussian states. This set will be abbreviated as LOG, and the 
statement that a transformation from a state p, pure or mixed, 
to a state p' is possible will be written as 

p — > p' under LOG. 

In quantum optical systems this class of operations is the one 
that can be realized (with present technology) as a combina- 
tion of applications of beam splitters, phase shifts, and squeez- 
ers together with the possibility to append additional field 
modes locally. 
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FIG. 1. Any local Gaussian completely positive map can be con- 
ceived as a composition of a local joint symplectic transformation on 
both the system and additional oscillators which have been prepared 
in a Gaussian state and a partial trace operation with respect to the 
additional oscillators. 

The physical system under consideration is a bi-partite 
quantum system with one canonical degree of freedom each, 
such as two modes of an electromagnetic field. As in Refs. 
[Q] such a system will be called a 1 x 1-system, consisting 
of parts 1 and 2. In order to exploit the elegant formalism 
that is available to describe Gaussian quantum states [||^] it is 
convenient to group the Hermitian operators corresponding to 
position and momentum in a vector, O = (Xi, Pi,X 2 , P^)- 
The canonical commutation relations (CCR) can then be sub- 
sumed into the skew symmetric block diagonal 4 x 4-matrix 
£ according to [O n ,O m ] = iT, nm , n,m — 1, ...,4. For a 
given state p, for which the second moments exist, let the real 
4 x 4-matrix T be defined as 
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Y nm = 2tr [p (0„ - (0„) P ) (O m - (O m ) P )] - iS 



flf(o, 6, c, d) := (a 2 - 1) + (6 2 - l)c 2 d 2 + 2cd - ab(c 2 + d 2 ). 



where (O n ) p — Xx[pO n ]. The matrix T will be referred to as 
covariance matrix. Not all symmetric 4 x 4-matrices are le- 
gitimate covariance matrices: the restriction that p is a state 
manifests itself as the condition r — iS > for the covari- 
ance matrix, which is in fact a formulation of the uncertainty 
relations. For Gaussian states Jl2[ ] the covariance matrix to- 
gether with the mean values of the position and momentum 
operators are sufficient to fully specify the state. The first mo- 
ments, however, are of no relevance for the issue of this paper, 
because they can always be made to vanish by an appropriate 
local translation in phase space. 

We will now turn to the possible state transformations. Of 
particular interest are the linear transformations from one set 
of canonical coordinates to another set which leave the CCR 
invariant. In a system with two canonical degrees of free- 
dom they form the group of real symplectic transformations 
Sp(4, R) [§]. The group Sp(A, R) consists of the real 4x4- 
matrices S obeying S T HS = S; the group Sp(2N, R) can 
be defined in an analogous manner for iV canonical degrees 
of freedom. Under a symplectic transformation a covariance 
matrix is transformed according to Y i — ► S T TS. On the 
level of states it is accompanied by a unitary operation p i — ► 
U (S)pU (Sy, then called symplectic operation. A local sym- 
plectic transformation is a matrix S of the form S = Si © S2, 
where Si, S2 £ Sp(2, R). The most general LOG can now 
be conveived as a composition of a joint symplectic transfor- 
mation S = Si © S 2 with Si, S 2 £ Sp(2N + 2, R) on the 
original systems 1 and 2 and on two additional systems with 
TV canonical degrees of freedom each of which has been lo- 
cally prepared in a Gaussian state, and a partial trace operation 
with respect to the additional systems (see Fig. 1). 

Any covariance matrix Y of a bi-partite 1 x 1-system can be 
written in block form as 

where Ai,A 2 , and B are real 2 x 2-matrices []13[], One can 
uniquely characterize the orbit 0(Y) of Y with respect to lo- 
cal symplectic transformations by a vector (6,6,6,6) £ 
R 4 , the entries of which are given by 6 := l^il 1 / 2 > 1, 
£2 := l^bl 1 / 2 > 1, where |.| denotes the determinant. £3 
and £4 are the solutions of 66 = \B\, £| + £f = (\B\ 2 — 
|T| + lAillAaD/dAxll^al) 1 / 2 , such that 6 > 161- It has 
been shown in Ref. [p[ that T can always be transformed into 
a covariance matrix S^YS which is of 'normal form' by using 
an appropriate local symplectic transformation S: this means 
that S T TS is of the form of Eq. (0), but with B = diag(6, 6) 
and Ai = diag(&,&), i = 1, 2. 

Whether a transformation of a state p to a state p' with 
respective covariance matrices Y and Y 1 is possible or not, 
will turn out to be largely determined by two functions 
fi^ T , f 2 ^ r that will be called minimal functions for rea- 
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sons that will become clear later. Let g : R + x E -* R, 



For a pair (Y, V) of covariance matrices with associated vec- 
tors (6 , £ 2 , 6 , 6 ) and j£ . f 2 » f 3 > Ci) with & , 6 > define 
the two functions /f , / 2 r_>r ' : R + x R — > R as 

A r ^ r := g(^,6,£/6,y/6), 

f 2 T ^ r '(x,y) :=g(e 2 , &,&/x,t,'Jy). 

The first statement concerns LOG in system 1 only. In this 
case the criterion amounts to simple inequalities that have to 
be satisfied. The second gives the full result for general LOG. 

Proposition 1. - Let p and p" be Gaussian states of a lxl- 
system with covariance matrices Y and Y" and associated 
vectors (£1,6, £3, £4) and (£', £2, £3, £4) with > 0. 

Then p — ► p" under LOG in system 1, if and only if 

1. 1661/6 > m'i\/tf, 2. /r r "(M) > 0. 

Proposition 2. - Let p and p' be Gaussian states of a lxl- 
system with covariance matrices Y and V and associated vec- 
tors (6 ,6,6,6) and (6 , £2 ? £3 > £4) w/f/1 £4,6 > °- r/ie " 
/? > p' under LOG, if and only if one of the points 

{x,y)e(fl^')-\Q)n{fl^')-\Q) 

satisfies |66l^/6 > > 

Proof of Proposition 1. We begin with investigating what con- 
ditions have to be met when a LOG is implemented in system 
1 and a symplectic operation in system 2. Starting point is 
a general representation theorem concerning Gaussian com- 
pletely positive maps [|TT|]: Any Gaussian completely positive 
map is reflected on the level of the covariance matrix as a map 

T 1 — ► M T YM + G, (2) 

where M and G are real 4 x 4-matrices, G is moreover sym- 
metric. The condition 

G + ffi - iM T HM > (3) 

on the matrices M and G incorporates the complete positivity 
of the map. The state transformation mapping Y on Y" can 
be decomposed into three steps: first, an appropriate matrix 
S = Si © S 2 , Si, S 2 £ Sp(2, R), is applied on the initial co- 
variance matrix V, such that S T TS is of normal form. Then a 
LOG restricted to system 1 and a symplectic operation in sys- 
tem 2 is implemented, mapping S T TS onto another matrix in 
normal form. Finally, T = T x © T 2 , T U T 2 £ Sp(2,R), is 
used in order to transform the resulting matrix into Y". The 
second step can be represented in the form of Eq. (Q) with 
real matrices M and G. Clearly, the composition of the three 
steps, T 1 — ► (T T M T S T )T{SMT) + T T GT amounts again 
to a LOG. Therefore, we can without loss of generality assume 
that both r and Y" are already in normal form with associated 
vectors (6,6,6,6) and (f» £» , &')■ The task is then 
to find appropriate real matrices M and G as above such that 
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r" = M T YM + G, representing a LOG restricted to system 1 
and a symplectic operation in system 2. Hence, it is required 
that M and G are of the form 



M = M 1 @M 2 ,G = G 1 ®Q, 



(4) 



where G\ is symmetric and Af 2 £ Sp(2, R), i.e., M 2 satis- 
fies M 2 T EAf 2 = S. Due to the normal form of V and T" 
we have that Af 2 r diag(^ 2 , £2)^2 = diag(£ 2 ,£2)i and it fol- 
lows that M 2 € 50(2). Let us set M 33 = M 44 = cos(0/2), 
M34 = -sin(0/2), and M43 = sin(0/2) with e (-2tt, 2tt]. 
The requirement that r" = M T TM+G implies then a certain 
set of equations that has to be satisfied, connecting the entries 
of Mi and M 2 . An elementary calculation yields finally 

Mn = (£ 3 7£ 3 )cos(0/2), M 22 - (£7£ 4 )cos(0/2), (5) 
M 12 = -(££/£ 3 )sin(0/2), Af 2 i = (&7&)sin(0/2). (6) 

Not all such matrices AT = Afi © A/ 2 and G = T" - M T YM 
define a completely positive map, however. Due to the block 
diagonal form of M and G, the inequality (j^) reflecting the 
complete positivity can be written as 

H 1 :=Gi-H£(l-|Mi|) >0. 

As Hi is a Hermitian 2 x 2-matrix, Hi > is in turn equiv- 
alent to l-ffil > 0, tr[i?i] > 0. The determinant \Mi\ = 
(£3 £4 )/(?3?4) i s independent of 0. The determinant and trace 
of Hi can be evaluated to tr[i?i] = 2f£' - £i||Afi|| 2 and 

\Hi\ = (e") 2 -ei£il|Mi|| 2 +£?|Mi| 2 -(l-|Mi|) 2 , where 
|| Mi || denotes the Hilbert-Schmidt norm of Mi, Hence, it is 
always optimal to chose in such a way that 



||Mi|| 2 = 



(£ 



ll\2 



£ 
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is minimal. But since £f > £| and (£3) > (£4 ) , it is true 
that always (£ 3 '/£ 3 ) 2 + ^{/U? < (&7£ 3 ) 2 + (£ 3 7^) 2 , and 
therefore, = is the optimal choice. To simplify the struc- 
ture of the requirements one can proceed as follows: The in- 
equality \Hi\ > implies in particular that 



£7€i > l&'&'l/lfck 



(7) 



Whenever (Q) is satisfied, |ifi| > yields a stronger upper 
bound for \\Mi || 2 as tr[Hi] > does, as then 



m) 2 +ex\Mi\ (1 - imid 2 )/(^i) < 2^76 



(8) 



Altogether, this implies that equivalently to requiring the va- 
lidity of both \Hi\ > and tr[i?i] > one may require that 
both Eq. (0) and \Hi\ > hold. Therefore, we finally arrive 
at the statement that p — > p" under under a LOG in system 
1 and a symplectic operation in system 2 if and only if both 
< 16641/6 and /f^ r > are satisfied. 
This criterion depends only on the invariants with respect to 
local symplectic operations, and hence, we arrive at Proposi- 



It is worth noting what physical situation is reflected by 
equality /f^ r "(£ 3 , £4) = 0. Under the constraint that the 
entries of M and G are given by Eqs. (f|,|j|6|), it can be shown 
easily that equality holds if and only if M and G satisfy 
G = K T G- Y K, where K := M T T,M - S. Solutions of 
this type are the minimal solutions in the sense of Ref . [Q] : 
For a given initial covariance matrix T and a given matrix M 
the equation G = K T G~ 1 K specifies those symmetric ma- 
trices G that add minimal noise. Such LOG will consequently 
be called minimal. Hence, a LOG in system 1 from p to p" is 
minimal if and only if ff ^ r (£3 , £4) = holds (see Fig. 2). 
Therefore, one may interpret the conditions of Proposition 1 
in physical terms as follows: the first condition requires that 
the 'stretching' \M\ is sufficiently small, the second makes 
sure that enough noise is introduced in the course of the trans- 
formation. 

Proof of Proposition 2. A general LOG can again be de- 
composed into several steps. As before, without loss of gen- 
erality one may assume that the initial and the final covariance 
matrices Y and Y' are of normal form with associated vectors 
(£1,6, £3, £4) and (£i, £2, £3, £4). respectively. In two inter- 
mediate steps one transforms Y 1— > Y" and Y" 1— > Y' by means 
of LOG restricted to one system and appropriate symplectic 
operations in the other system. The vector associated with V" 
will be denoted as £2, x, y). One can proceed as before, 
and after applying analogous steps one finally arrives at the 
criterion that p — ► p' under LOG if and only if there exists 
an (x, y) £ R + x R such that the inequalities 



wi > \xy\ > w 2 , fi^ r '(x,y) > 0, /. 



2 ^ r '(x,y)>0 (9) 



tion 1. 



□ 



are simultaneously satisfied, where Wi := |£ 3 £4|£J/£i and 
w 2 : = l£ 3 £4l£2/£ 2 - This is already a criterion of its own, 
but it still requires a search in a two-dimensional set. The 
key observation in a simplification is that the intersection of 
the interior of the set (f[^ r )~ 1 (0) and and the set Ni := 
{{x,y) G R + x R||a;?/| = w{\ is empty for both i = 1,2 
(see Fig. 2). The minimal value of (cc/£ 3 ) 2 + (y/£4) 2 for 
(x, y) S Ni is given by 2£(/£i. Hence, it follows from Eq. 
(§) that fi(x,y) < for all (x,y) € N x , if £1 ^ and 
fi(x,y) < for all (x,y) € N x , if £1 = ££. Similarly, 
f 2 (x,y) < for all (x,y) £ N 2 . Moreover, /f^ r is con- 
tinuous on R x R+, and for f 2 ^ r there exists a continuous 
continuation on R + x R. The problem is therefore reduced 
to the subsequent search for intersection points: there exists 
an (x, y) G R + x R satisfying if and only if there ex- 
ists a point (x,y) e (fi^ r 'y 1 (0)r\(f^-* r 'y 1 (0) such that 
wi > \xy\ > w 2 . This is Proposition 2. In particular, this 
means that if the transformation p — ► p' is possible, it can 
always realized as a composition of two minimal LOG in sys- 
tem 1 and 2, respectively Jiff], rj 
So far, the simple case has been omitted that the initial 
state p has a covariance matrix Y with associated vector 
(£i) £2, £3, £4), where £ 4 = 0. It turns out that one can pro- 
ceed as before. In the notation of Proposition 1 (but with 
£4 = 0), one arrives at the statement that p — ► p" under 
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LOG restricted to system 1 , if and only if both £" = and 
(#/6) 2 < m) 2 ~ Consequently, in the no- 

tation of Proposition 2, p — > p' under LOG, if and only if 

i'l = and (£'/6) 2 < ((^) 2 - 1)((^) 2 - 




-1 1 
FIG. 2. Given is a state p with a covariance matrix Y with associ- 
ated vector (£i, £2, £3, £4) = (3, 5, 1, 1/2). The shaded area depicts 
what values of £3 and £4 are accessible under a LOG in system 1, un- 
der the assumption that the final covariance matrix V" is associated 
with a vector (£1,62, £3 ,64) with £" = 2. The thick line corre- 
sponds to those points (£3 , £4) with /f^ r " (£3, £4) = for which 
the transformation is a minimal LOG, the dashed line represents the 
points satisfying = |£ 3 £4|£i7£i- 

As a first application we can look for Gaussian states p and 
p' that are incommensurate, that is, pairs of states (p, p') for 
which neither p — > p' under LOG nor p' — > p under LOG 
holds. As can readily be verified using Proposition 2, an ex- 
ample of such a pair is given by states specified by covariance 
matrices associated with (2, 2, 1, 1), (2, 2, 1, —1/2), respec- 
tively. The relation that a state can be transformed into an- 
other state under LOG induces hence a partial order on the set 
of Gaussian states, but not a total order. 

With this letter we have posed and answered a basic ques- 
tion: Under the constraint of locality, we ask which pairs of 
Gaussian states allow for a transformation from one state to 
the other. The choice for the set of allowed operations - Gaus- 
sian completely positive maps - has been motivated by prag- 
matic considerations: in quantum optical systems such opera- 
tions can be implemented with present technology. Needless 
to say, there are many open questions that may be approached 
with similar methods. In particular, one may take into ac- 
count selective measurements projecting on Gaussian states 
together with classical communication. It is the hope that this 
letter stimulates such further considerations. 
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